Ideal glass-glass transitions and logarithmic decay of correlations in a simple system. 
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We calculate the ideal-glass-transition line for adhesive hard spheres in the temperature- volume- 
fraction plane within the framework of the mode-coupling theory. We find two intersecting lines, 
controlled by the hard-core and the adhesive part of the potential respectively, giving rise to two 
different mechanisms for structural arrest. In the glass region we identify the presence of a glass- 
glass-transition line ending in a cusp bifurcation which causes, even in the close by liquid region, a 
logarithmic decay of correlations. 
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The crossover from a liquid to an amorphous solid, 
observed near the calorimetric glass transition tempera- 
ture Tg, exhibits as a precursor phenomenon an anoma- 
lous dynamics, called glassy dynamics. Its evolution is 
connected with a critical temperature Tc above Tg. It 
has been studied extensively in the recent literature of 
the glass-transition problem, both experimentally [Q-^, 
numerically ||^J^ and theoretically Experiments 
around Tc have been interpreted in the frame of the 
mode-coupling theory (MCT) for structural relaxation. 
MCT deals primarily with closed equations of motion for 
the normalized density-fluctuation-correlation functions 
^q{t) for wavevector moduli q. The equilibrium structure 
enters as input in these equations via the static structure 
factor Sq. The theory explains Tc as a glass-transition 
singularity resulting as a bifurcation phenomenon for the 
self-trapping problem of density fluctuations. Below Tc 
the interaction of density fluctuations leads to arrest in 
a disordered array, characterized by a Debye- Waller fac- 
tor fq > 0. Near the transition, the MCT equations can 
be solved by asymptotic expansions. Beyond the initial 
transient dynamics, correlation functions are predicted to 
decay with a time fractal toward a plateau value f^, the 
critical Debye- Waller factor. Above Tc the correlations 
decay from to zero, and this is the MCT interpreta- 
tion of the a-process of the classical literature of glassy 
dynamics jl^. The initial part of this decay is another 
time fractal, called von Schweidler's law. The values of 
the fractal exponents are controlled by the so-called ex- 
ponent parameter A < 1, which depends solely on Sq. For 
details and citations of the original literature the reader 
is referred to the review in Ref. A number of tests 
of MCT results against data, among them the ones in 
Refs. p]-[7|, pT| , p^ , demonstrates that this theory treats 
reasonably the evolution of structural relaxation in some 



systems. 

The MCT bifurcations are caused by a nondegenerate 
eigenvalue of a certain stability matrix to approach unity 
from below . Therefore the bifurcation scenario for fq 
is that known for the zeroes of a polynomial as induced by 
changes of the polynomial's coefficients [Q . The generic 
case for a change due to a single control parameter is a 
fold bifurcation i.e. the coalescence of two real ze- 
roes to a degenerate one. In all the microscopic models 
for liquid-glass transitions studied so far, only the fold bi- 
furcation singularity had been identified. The next more 
complicated scenario is the cusp bifurcation, equivalent 
to the coalescence of three real zeroes. Generically, this 
case can be obtained only if two control parameters are 
varied . The MCT-bifurcation dynamics for this case 
has been worked out in a leading order asymptotic expan- 
sion . It is drastically different from the one for a fold 
bifurcation in the sense that relaxation stretching is much 
more pronounced. In particular there appear logarithmic 
decay laws, i.e. l//-noise spectra. Cusp bifurcations in 
MCT are endpoints of transition lines characterized by 
A = 1, and typically they are located near intersections 
of bifurcation hues Sjogren identified some sets 
of polymer dielectric-loss spectra which could not be fit- 
ted by the fold-scenario but were compatible with a cusp 
bifurcation. However, in the data used in Ref. jl6j the 
temperature was varied as the only control parameter, 
and therefore a complete test of the cusp-bifurcation sce- 
nario was not possible. 

In this paper we want to demonstrate that a cusp bifur- 
cation near the intersection of two fold-bifurcation lines 
is possible in a simple but realistic model, if the two con- 
ventional control parameters temperature T and density 
n are varied. With this aim we study an extension of the 
hard sphere system (HSS), the sticky hard spheres sys- 
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tern (SHSS) introduced by Baxter ||T^ and used exten- 
sively, in particular in colloid physics The HSS 
is the archetype of a simple system and its only control 
parameter is the packing fraction = 7rncr^/6, with a 
denoting the hard sphere diameter. It has been studied 
by dynamic light scattering in its realization as a certain 
colloidal suspension. The system exhibits an ideal glass 
transition at some critical packing fraction (^^ |4) . This is 
a sharp transition of the fully equilibrated sample from 
an ergodic liquid to a nonergodic amorphous solid. The 
glassy dynamics of the colloid exhibits the same scaling 
laws and fg anomaly as found for the more conventional 
systems |jlp^,^,0,0|l^] . The measured correlators can 
be explained quantitatively by the first-principle calcu- 
lation produced by MCT [||. The SHSS potential is the 
limiting form of a square well, infinitely narrow and in- 
finitely deep ll^. It is defined by the interparticle po- 
tential i5V{r) = oo for r < cr, ln{12TA/{a + A)) for 
a < r < a + A, and for r > tr -I- A, in the hmit A ^ 0. 
Here (3 — {kBT)~^, ks is Boltzmann's constant, and r 
is the adhesiveness parameter which plays the role of a 
dimensionless temperature. For t ^ oo one recovers the 
HSS. The Ornstein-Zernike equation Sq = 1/(1 — ncq) 
with the Percus-Yevick closure has been solved analyti- 
cally by Baxter and we use his expressions for Sq, as a 
function of r and (f) JTzt as input of this work. The cru- 
cial point is that the interaction exhibits two features, 
which favour arrest of density fluctuations for two quite 
different reasons. One is a strong short-ranged repulsion. 
It is the mechanism leading to the cage effect for liquid 
dynamics and thereby to arrest driven by density fluctu- 
ations near the structure factor peak position. The other 
one is a short-ranged attraction. It leads via the usual 
mean-field mechanism to a softening of the density fluc- 
tuations, in particular of those for small wave vectors. It 
is the mechanism which leads to a liquid-gas transition 
at low values of r. 

The MCT equations of motions are j|] 



^q{t)+Vq^q{t)+VLi^^q{t)+ni^ / ds TO, - s) (s) = 



(1) 



Here Vlq — ^Jqv j Sq, with v denoting the thermal velocity, 
is an effective phonon-dispersion law and Vq ~ viq^ de- 
notes a damping constant (we use vi = 1). The kernel niq 
is given as mq{t) — !Fq[<^k(t)], where the mode-coupling 
functional Tq is determined by the structure factor: 



[h] - \ 



(2^ 



q-kck + q- [q- fc)c|g._fc| 



(2) 



(3) 



spaced q values extending up to qa = 72 fq is ob- 
tained by an iterative solution of the bifurcation equation 



fq/il- !q)=Tq[h]. 



(4) 



We numerically solved Eq. on a grid of 400 equally 



For T oo we recover the known HSS result, for which 
the ideal glass transition at (f)^^^ ^ 0.516 [|| is lead by 
the increase of —Cq and Sq with (f). For > (j>^^^ , the 
particles are trapped in cages formed by their neighbors. 
On decreasing r attractive forces become relevant; near- 
est neighbor pairs are closer and therefore holes in the 
cage are produced, destabihzing the glass. A decrease of 
r has to be compensated by an increase of (j) and thus the 
Tc versus (jjc high-temperature-transition line has to bend 
away from the (j) = <j)^^^ value. This effect is brought out 
by the MCT solutions as shown by the dashed lines Bi in 
Fig. |l|. The destabilization of the glass on decreasing r is 
related, within MCT, to the increase of pair correlations 
which weaken the screened potential —Cq and decrease 
the value of Sq^^^ , as shown in Fig. ||-top. Moving along 
El in Fig. [l] by decreasing r implies an increase of the 
attraction and thus of the particle localization compared 
to the HSS case. Consequently, increases and spreads 
out further in the g-domain as demonstrated in Fig. ^ 
top. 

The increase of Cq and Sq at small q becomes in- 
creasingly relevant on decreasing r producing a differ- 
ent mechanism for structural arrest. This increase of the 
compressibility is a precursor of the liquid-gas critical 
point. The enhanced amplitudes for density fluctuations 
increase their anharmonic interactions, and indeed the 
solution of the MCT equations predict their spontaneous 
arrest on a low-temperature-transition line B2 in the r vs 
(f> plane, shown in Fig. ^ Increase of suppresses small 
q fluctuations. But it also increases —Cq and the first 
peak of Sq , as shown by a comparison of results shown in 
Fig. ||-middle. Both effects stabilize the glass, since the 
mode-coupling constants, Eq. (|^) increase. This explains 
the upward bending of the full line in Fig. ||. Since ar- 
rest for small q is favored on the B2 branch, is much 
larger there than for the HSS. The decrease of the com- 
pressibility upon increasing (j) make spontaneous arrest 
less easy and therefore fq decreases upon compressing. 
Both features are demonstrated in Fig.|^middle. 

The explained interplay of attraction, and repulsion 
places the meeting point C* of the two branches Bi and 
B2 to a packing fraction 0* = 0.5565, higher than cj)^^^ , 
and to T* = 1.320. Thus there is a reentry phenomenon 
if T is lowered, e.g. for a packing fraction between 0.52 
and 0.55. At some r above 1.32 the glass melts since the 
increase of attraction-induced-pair-correlations destabi- 
lizes the cages for particle arrest. At some lower r the 
liquid freezes again into a glass, because of the compress- 
ibility increase. At C* , the fq calculated on the two 

different branches, /g^'' and are different, while the 
corresponding Sq are indistinguishable. Thus two types 
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of glasses, differing in their fq and in their dynamical 
properties, can be obtained close to C* . In the region 
where both /g^^ and fq"^"^ are solutions of the Eq. (jj), 
the one with the larger fq is the relevant solution. This 
exemplifies the maximum theorem of MCT: if there is 
a solution of Eq. (Q), say fq, the long time limit fq at 
the same equilibrium state obeys fq > fq For our 

c(2) 

case it means that the continuation of fq takes over 
the role of fq, i.e. Bi stops at C*. Therefore the two 
branches of the transition line do not join smoothly at 
C* . The i?2 branch continues, as shown by the inset of 
Fig. until it reaches a cusp bifurcation as an endpoint 
of the transition line. Indeed, we evaluate the exponent 
parameter A along the Bi and B2 lines and find, as shown 
in the inset of Fig. [^bottom, that along the B2 line A 
keeps on increasing till it reaches unity for r = 1.37 and 
(j) = 0.55807. The line between C* and the cusp point is 
a line of glass-to-glass transitions. For the given equilib- 
rium structure for (r, (jj) on this line, there are two pos- 
sibilities for structural arrest. The one on high r side is 
driven by the excluded volume mechanism and the other 
one by the attraction-induced compressibility increase. 
The latter mechanism leads to a larger fq than the for- 
mer. The endpoint is characterized by their difference 
approaching zero, as shown in Fig. [^bottom. 

In Fig. ^ the glassy dynamics outside the transient is 
exhibited for states in the reentry region. The dynamics 
is, up to a regular time scale, independent of the tran- 
sient. In particular it is the same for a colloid, as for a 
conventional liquid for which Eq. (|l|) is formulated [§. 
The uppermost curve corresponds to a point in the glass 
region and demonstrates arrest near the Bi branch, the 
others exhibit stretched liquid relaxation to zero. In- 
terestingly enough, the exhibited anomalous dynamics 
is dominated by the cusp singularity, not by the closer 
fold. Indeed even though a huge dynamical window is 
considered, the known fold bifurcation pattern cannot 
be recognized. There is no power-law decay towards the 
fq^^^ nor is there an a-process obeying the superposi- 
tion principle. We checked that the mentioned asymp- 
totic fold-bifurcation features appear only after tuning 
much closer to the transition line, which implies the ex- 
tension of the dynamical window to even larger sizes. 
Instead, the cusp dynamics appears in the liquid re- 
gion. The stretched relaxation in the three-decade win- 
dow < logio{t /to) < 3 shows the approach towards the 
critical decay of the cusp . Then the correlators fol- 
low closely a logarithmic law in time. For = 0.555 it 
extends over 5 decades. Upon decreasing (j) its range of 
validity shrinks but even at </> = 0.540 it extends over two 
decades. This logarithmic decay is the known cusp sub- 
stitute for the start of a fold a-process. Let us emphasize 
that the above scenario does not require a fine tuning of 
T. For r = 1.32 the pattern is similar to the one shown 
for T = 1.40. 



In summary, we have shown that the SHSS presents 
a peculiar structural arrest dynamics at high packing 
fractions due to the competing mechanism of hard core 
and attractive interactions. Two differently sloped ideal- 
glass transition lines appear in the phase diagram. In 
the region where these two lines meet the slow dynam- 
ics changes from power-law to a logarithmic-law due to 
the influence of a nearby cusp singularity. In most tests 
of MCT the data had been compared with the univer- 
sal results obtained for the asymptotic dynamics near a 
fold transition. But our results show that even a sim- 
ple system can exhibit glassy dynamics, which does not 
exhibit asymptotic laws within accessible dynamical win- 
dows. It seems relevant to test, by spectroscopy for col- 
loidal systems and by molecular-dynamics simulations, 
whether specific systems exhibit the shown complex dy- 
namics and, in case, whether it is handled properly by 
the present MCT calculation. 
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FIG. 1. Ideal glass-transition lines for the Baxter model 
(squares and dots). Bi (dashed) is the high temperature 
branch of the glass-transition line and B2 (full) is the low 
temperature one. The arrow indicates the critical packing 
fraction of the hard sphere system. The open circle denotes 
the cusp-point. 
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FIG. 2. Structure factors Sq/10, direct correlation func- 
tions Cq and critical Debye- Waller factors along Si (top), 
B2 (middle) and along the glass-glass transition lino (bottom). 
Top: Tc = 8 , (j)c ~ 0.530 (light dashed lines) and Tc — 1.32, 
(p^ = 0.522 (heavy dashed lines). The dotted line in the in- 
set is Cq at 0f = 0.516 for r = oo. Middle: Tc = 0.2, 
0c = 0.295 (light fuU lines) and Tc = 1.31, 0c = 0.5561 (heavy 
full lines). Bottom: The full (dashed) lines refer from top to 
bottom (bottom to top) to t = 1.34, 1.35, 1.36 respectively, 
on the B2 {Bi) side of the glass-glass transition line. The 
heavy full line is fq for the cusp point t = 1.37, = 0.55807. 
The lower Sq refers to the same states. The inset shows the 
exponent parameter A along Bi (dashed) and B2 (full). 
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FIG. 3. Density correlator for qa = 14.4 at r = 1.40 and 
various as a function of logio{t/to) with to = cr/v. 
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